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ABSTRACT 

We identify a very deep sample of 198 extremely red objects (EROs) in the Chandra 
Deep Field South, selected on the basis of 1775 — Ks > 3.92, to a limit Kg — 22 using 
the public ESO/GOODS survey. 

The ERO number counts flatten from a slope of 7 ~ 0.59 to 7 ~ 0.16 at K > 19.5, 
where they remain below the predictions for pure luminosity evolution, and fall below 
even a non-evolving model. This suggests there is a significant decrease with redshift 
in the comoving number density of passive/very red galaxies. 

We investigate the angular correlation function, uj{9), of these EROs and detect 
positive clustering for Kg = 20.5-22.0 sources. The EROs show stronger clustering 
than other galaxies at the same magnitudes. The uj{9) amplitudes are best-fitted by 
models in which the EROs have a comoving correlation radius tq ~ 12.5 ± 1.2 
Mpc, or ro — 21.4 ± 2.0 Mpc in a stable clustering model. 

We find a 40-arcsec diameter overdensity of 10 EROs, centered on the Chandra 
X-ray source (and ERO) XID:58. On the basis of colours we estimate that about 7, 
including XID:58, belong to a cluster of EROs at z ~ 1.5. 

The 942 ksec Chandra survey detected 73 X-ray sources in the area of our ERO 
sample, 17 of which coincide with EROs. Of these sources, 13 have X-ray properties 
indicative of obscured AGN, while the faintest 4 may be starbursts. In addition, we 
find evidence that Chandra sources and EROs are positively cross-correlated at non- 
zero (~ 2-20 arcsec) separations, implying that they tend to trace the same large-scale 
structures. 

In conclusion these findings appear consistent with a scenario where EROs are 
the z > I progenitors of elliptical/SO galaxies, some forming very early as massive 
spheroids, which are strongly clustered and may evolve via an AGN phase, others 
more recently from mergers of disk galaxies. 

Key words: galaxies: evolution; galaxies: elliptical and lenticular, cD; galaxies: high- 
redshift; X-ray: galaxies 



1 INTRODUCTION 

The 'extremely red objects' (EROs) are a population of very 
red galaxies (e.g. with R — K > 5 ot I — K > 4) which 
appear faintward of /C" ~ 17.5. EROs are of special inter- 
est in the study of galaxy evolution, in that their colours 
and other properties suggest that they are the high-redshift 
(2 ~ 1-2) counterparts and progenitors of local elliptical and 
SO galaxies, and are amongst the oldest galaxies present at 
these redshifts. 

Recent spectroscopy of some of the brightest EROs 
(Cimatti et al. 2002a; Saracco et al. 2003) has revealed 
a mixture of old, passive galaxies (hereafter pEROs) and 
younger star-forming or post-starburst galaxies with strong 



dusty-reddening (hereafter dsfEROs). Some EROs also host 
AGN activity. Recent Chandra surveys find that ~ 10-30 
per cent of EROs on the Chandra Deep Field North (CDFN) 
are detected in X-rays and, while in some of these the X- 
rays may be produced by starbursts, the majority have X- 
ray properties indicative of obscured AGN (Alexander et al. 
2002; Vignali et al. 2003). 

In our previous paper (Roche et al. 2002, hereafter 
also Paper I), we combined our deep K-band survey of the 
'ELAIS:N2' field with the UFTI (reaching K ~ 21) and In- 
grid {K ~ 20) cameras with a deep 7?-band image, to select 
a sample of 158 EROs as R — K > 5.0 galaxies. The num- 
ber counts of these EROs were lower than predicted by a 
model in which all E/SO galaxies evolve by 'pure luminosity 
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evolution' (PLE), but could be fitted by a model combin- 
ing passive luminosity evolution with galaxy merging and a 
decrease with rcdshift in the comoving number density of 
passive galaxies ('merging and negative density evolution'). 

A number of studies have found that EROs are more 
strongly clustered than other galaxies of similar magnitudes 
(Daddi et al. 2000; Firth et al. 2002). As a measure of cluster- 
ing, we investigated the angular correlation function, u!{9), 
of our ELAIS:N2 sample. Our results, combined with those 
of Daddi et al. (2000) and Firth et al. (2002), were best-fitted 
by models in which the EROs have a comoving correlation 
radius ro ^ 10-13 h~^ Mpc - which implies even stronger in- 
trinsic clustering than that of present-day L > L* ellipticals 
(ro ~ 8h-^ Mpc). 

We examined the morphology of the brighter {K < 
19.5) EROs (31) and estimated that ~ 60 per cent were 
bulge-type (E/SO) galaxies with the others a mixture of 
disk, irregular and merging galaxies. This mix is consistent 
with the WFPC2 study of Moriondo, Cimatti and Daddi 
(2000). The angular sizes of the EROs were generally con- 
sistent with our best-fitting evolution model, whereas a non- 
evolving model underpredicted their surface brightness. 

Seven of these 31 brighter EROs were detected as radio 
sources on a deep VLA survey to ^(1.4 GHz) ~ 27.6/iJy. 
Six were point-like within resolution 1.4 arcsec. One of these 
was a powerful radio galaxy, and the other five were much 
fainter radio sources within elliptical galaxies, which could 
have been either obscured central starbursts or obscured 
AGN. 

The seventh was an elongated radio source, aligned 
with an extended, low surface brightness optical morphol- 
ogy, suggestive of a dusty starburst galaxy. This implied that 
some (but perhaps only a few per cent) of EROs have star- 
formation rates as high as the > 100-1000 Moyr~^ needed 
for radio detection in this survey. Small et al. (2002) detected 
11 EROs on an even deeper VLA image, and found several 
other examples of dusty starburst galaxies with extended 
and aligned radio emission. 

A scenario was proposed in which the (i) pEROs are 
strongly clustered primordial ellipticals which had formed as 
powerful sub-nini galaxies at even higher redshifts, (ii) the 
younger dsfEROs are starburst and recent post-starburst 
galaxies formed from the merging of disk galaxies, which 
would only later become passive E/SOs. At all epochs, inter- 
actions cause new dsfEROs to be added to the ERO class, 
increasing the comoving number density (as in our M-DE 
model). 

If the progenitors of dsfEROs are less clustered than the 

old pEROs, this process would be accompanied by a dilution 
of ERO clustering with time (i.e. an increase with redshift). 
Daddi et al. (2002) have recently claimed that spectroscop- 
ically selected pEROs axe indeed more clustered than ds- 
fEROs. However, their sample was small (33) and included 
only the brightest EROs (mostly K < 19.2), and the clus- 
tering was analysed in redshift space only, so this needs to 
be verified with larger and deeper samples. 

Another result of interest is the Almaini et al. (2003) 
finding that sub-mm (detected with SCUBA) and X-ray 
(Chandra) sources on the ELAIS:N2 field, were positively 
cross-correlated at non-zero separations. This was inter- 
preted as possible evidence for an evolutionary connection 
between the two. 



In this scenario (Archibald et al. 2002; Granato et al. 
2003), the sub-mm sources are giant ellipticals in the pro- 
cess of formation. At early epochs, the most massive galaxies 
would form stars very rapidly, producing far-infra-red flux 
(redsliifted to the sub-mm). In parallel with this, a super- 
massive, dust-enshrouded black hole is formed in the centre 
of the galaxy. After ~ 1 Gyr, the emission from the accreting 
black hole becomes sufficient to eject the dust and remain- 
ing gas from the galaxy, halting further star formation, and 
the galaxy ceases to be sub-mm bright but instead is trans- 
formed into an X-ray luminous QSO. 

The relevance of this to our present study of EROs is 
that, at the end of the QSO phase, the massive, gas-stripped 
galaxies can only become passively evolving red ellipticals, 
i.e., EROs. One consequence of this might be a similar cross- 
correlation between EROs and unreddened, non-ERO Chan- 
dra sources. 

The recent public release of the first installment of 
the ESO/GOODS survey - with optical and near- infra- red 
imaging of the Chandra Deep Field South (CDFS) pro- 
vided us with an excellent opportunity to investigate both 
the counts and clustering of EROs (to a fainter limit than 
previously) and their association with a very deep sample of 
X-ray sources. 

The layout of this paper is as follows: Section 2 describes 
the observational data and any data reduction. Section 3 the 
identification and colour-classification of the EROs, Section 
4 the ERO number counts, comparing those with evolution- 
ary models. Section 5 the ERO angular correlation function, 
Cij{0), with an interpretation in terms of intrinsic cluster- 
ing. In Section 6 we investigate the numbers of close pairs 
of EROs and identify any larger clusters, and in Section 7 
identify EROs which are also Chandra X-ray sources and 
look for a correlation between these two classes. Section 8 
concludes the paper with a discussion of our findings. 

In Paper I we had assumed a cosmology with Qm ~ 0.3, 
f^A = 0.7 (modelled using the analytic form of Pen 1999) 
and Ho = 55 km s~^Mpc~^. In view of new results from 
the Wilkinson Microwave Anisotropy Probe (e.g. Spergel et 
al. 2003), we here retain this Qm/^k, but adopt the higher 
Hi) — 70 km s~^Mpc~^, which reduces the time since the Big 
Bang from 17.16 to 13.48 Gyr. Some Ho-dependent quanti- 
ties are given in units oi h = {Ho)/WO km s~^Mpc~^ ) 



2 OBSERVATIONAL DATA 

2.1 Near Infra-red (JHKs) Imaging 

The Great Observatories Origins Deep Survey (GOODS) 
(Dickinson and Giavalisco 2003) is a public, multiwavolength 
survey of the Hubble Deep Field North and the CDFS. The 
ESO contribution to GOODS will include deep near-IR and 
optical imaging of the whole 150 arcmin^ of the CDFS. The 
noar-IR observations are being carried out using the 'ANTU' 
VLT (Very Large Telescope) , with the Infrared Spectrometer 
and Array Camera (ISAAC). ANTU was the first of four 8m 
telescopes installed at the European Southern Observatory 
(ESO) at Cerro Paranal, Chile. ISAAC is equipped with a 
Rowell Hawaii 1024 x 1024 pixel Hg:Cd:Te array covering 
l-2.5/L(?n, with a pixelsize 0.1484 arcsec giving a 2.5 x 2.5 
arcminute field-of-view. 
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ISAAC observations will eventually cover the whole 
CDFS, as a 32-frame mosaic in the J {Xcentrai = 1.25/xm) 
H (1.65/^m) and Ks (2.16/xm) passbands. ISAAC observa- 
tions of the first 8 GOODS fields, covering a contiguous area 
near the centre of the CDFS (approx. RA 3''32'"30^ Doc. 
—27 : 47 : 30) were completed on 1 February 2002 and made 
publicly available on 9 April 2002. 

2.2 7776-band Imaging with the HST ACS 

Another component of the GOODS survey is deep optical 
imaging of the CDFS with the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST). The ACS pro- 
duces a 4096 X 4096 pixel image with pixel size 0.05 arcsec. 
We make use here of newly released data from four phases 
of ACS observation, obtained from 3 October to 22 Decem- 
ber 2002, which cover the CDFS in mosaics of either 15 or 
16 fields. The CDFS imaging was performed in four ACS 
passbands {BVIZ). In this paper we make use of only the 
/77B-band imaging {\pivot = 0.769/im), where the exposure 
time per field is 1040 seconds for each phase. 

2.3 Chandra X-ray Observations 

We also make use of the CDFS X-ray source catalog of Giac- 
coni et al. (2002). The CDFS was observed with the Chandra 
X-ray Observatory for eleven pointings during 1999-2000, 
which were combined to give 942 ksec exposure time. The 
Chandra aim-point of RA 3''32'"28^ Dec. -27 : 48 : 30, 
where sensitivity and resolution are optimal, is approxi- 
mately concentric with our ISAAC data. 

Giacconi et al. (2002) detected a total of 346 sources 
in a soft (0.5-2.0 keV) and/or in a hard (2-10 keV) X-ray 
band, reaching aim-point detection limits -F(0. 5-2.0 keV) = 
5.5 X 10"" and F{2-1Q keV) = 4.5 x lO'^'' erg s'^cm'^, 
with resolution FWHM ~ 0.7 arcsec. Of these sources, 73 
are within the area of our ISAAC data. 



3 REDUCTION AND ANALYSIS 
3.1 ISAAC JHKs Data 

The 8 fields of ISAAC data were obtained from ESO 
(www.eso.org) in a reduced form, as 24 images (8 fields x 3 
passbands) each with a weighting map, astrometry and pho- 
tometric zero-point. Total exposure times in each field were 
20.8-28.6 ksec in K^, 13.92-18.12 ksec in H and 7.92-15.12 
ksec in J. 

For our purposes the only further processing required 
was to mosaic the 8 fields into a single image for each of the 
three passbands. However, there was considerable variation 
between these in signal-to-noise and in zero-point, and to 
take this into account: 

(i) we scaled the images and their weighting maps to a 
common photometric zero-point, that of the first field (the 
scaling of a weighting map is the inverse of that applied 
to the data image). Zero-points were also corrected for the 
Galactic extinction on the CDFS, but this is only 0.008 (J), 
0.005 {H) and 0.003 [K,) mag. 

(ii) We multiplied the images by their weighting maps. 
The resulting products were then mosaiced using iraf 



'combine' with 'offsets=wcs' (using the supplied astrome- 
try which is sufficiently accurate), and 'scale=none'. Note 
that this averages input images where they overlap in axea. 

(iii) The weighting maps were mosaiced in the same 

way. 

(iv) the mosaic from (ii) was then divided by the mosaic 
from (iii). The result of this is a 'final' mosaic where the 
contribution at each pixel from each of the input images is 
weighted by its appropriate weighting map. The division by 
the mosaiced weighting map normalizes the sum of these 
input image weightings to unity at each pixel. Hence the 
zero-point of the final mosaic should be uniform across all 
pixels and the same as that of the input images. 

(v) We still need a weighting rnap for the final mosaic. 
To obtain this the mosaiced weighting map from (iii) , repre- 
senting the average depth of the input images contributing 
at each pixel, is multiplied by a map of the number of these 
contributing images, thus giving a map of the total depth of 
observation at each pixel. 

The final mosaics for each of the three passbands, with 
their respective weighting maps, were then ready to be used 
for source detection. The mosaicing procedure weis checked 
by measuring source magnitudes on individual unmosaiced 
frames and comparing with those from the mosaic, and by 
checking for any sign of positional mismatch in overlap re- 
gions. 

3.2 ACS /775-band Data 

The ACS data was obtained from the HST archive 
(archive. stsci.cdu/hst/goods. html), already flat-fielded and 
calibrated to a common zero-point. The first and third (and 
the subsequently released fifth) observation phase covered 
the CDFS in a 3 X 5 mosaic of 15 pointings, with the long 
axis on position angle —22.35°. However, the second and 
third phase covered the CDFS in a 16 field grid orientated 
at 45° to this, ie. with position angle —67.35°. 

The number of pixels was too great for the whole ob- 
served area to be mosaiced into a single image. Instead, the 
four phases of observations were combined into 15 frames, 
corresponding to the 15 field areas of the first phase of ob- 
servation. Because of the two different orientations of the 
data, for each of these 15 frames there were typically 8 ob- 
served images which covered all or part of the field area. 
These data were combined by much the same methods as 
the ISAAC data above, but with an additional step of re- 
binning: 

(i) For each of the 15 frames, iraf 'wcsmap' was used to 
fit transforms between the pixel grid of the phase one obser- 
vation and that of all subsequent observations (typically 7) 
which cover at least part of the same area of sky. IRAF 'geo- 
tran' could then be used to rebin these observations into the 
same pixel grid as the phase one observation (which was not 
rebinncd). This rcbinning would allow observations taken at 
different position angles to be simply added together. 

(ii) The exposure maps supplied with each observed im- 
age were similarly rebinned into the phase 1 pixel grid. 

(iii) For each of the 15 frames, the data images to be 
combined (rebinncd, except for phase 1) were multiplied by 
their respective exposure maps and then added together us- 
ing IRAF 'combine'. The exposure maps were similarly com- 
bined. Then the combined datax exposure map image was 
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divided by the combined exposure map, to produce a 'fi- 
nal' co-added image where the contribution at each pixel 
from each of the input images is proportional to its weight- 
ing (from its exposure map), and the sum of the weighting 
factors is normalized to unity. 

(iv) A weighting map for each of the 15 'final' frames 
was produced by multiplying the combined weighting map 
by a map of the number of input images contributing at each 
pixel (again, as for the ISAAC data). 



3.3 Source Detection 

Sources were detected on the ISAAC and ACS images us- 
ing 'SExtractor'(Bertin and Arnauts 1996), as in Paper I. 
We make use of both 'total' magnitudes, derived by 'SEx- 
tractor' by fitting elliptical apertures to each detection, and 
aperture magnitudes, measured by 'SExtractor' in circu- 
lar apertures of fixed 2.0 arcsec diameter. The former will 
be used for limiting the sample and the latter for measur- 
ing colours. Magnitudes arc given in the Vega system and 
can be converted to the AB system as (Ks,H, J,I775)ab = 
iK„H, J, l775)vega + (1-841, 1.373, 0.904, 0.403). 

Source detection was performed on the ISAAC Ks mo- 
saic (together with its weighting map), with the chosen cri- 
terion that a source must exceed lAasky above the back- 
ground, or 22.75 Kg mag arcsec"^, in at least 6 contiguous 
pixels. In addition, a detection filter (3.0 pixel FWHM Gaus- 
sian) was employed. Detections on the low signal-to-noise 
edges of the mosaic wore excluded, leaving a data area of 
50.4 arcmin^. The catalog of Kg detections forms the source 
list for our ERG sample. 

Photometry of these sources in H and J was obtained 
simply by positionally registering (to the nearest pixel) the 
H and J mosaics to the Ks mosaic and then running SEx- 
tractor in 'double-image mode' to detect sources in K^, as 
before, and then measure their fluxes on the H and J images. 
Sources were detected separately on the ACS /rrs-band im- 
age with the criterion that they exceed 1.75asky in at least 
8 pixels, and a 2.5 pixel FWHM Gaussian detection fllter. 



3.4 Ks-baxid Galaxy Counts 

As some detections, especially at brighter magnitudes, will 
be Galactic stars, we performed star-galaxy separation on 
the source list, using a plot of peak/total Ks flux against 
magnitude. This was efTective to a limit Kg = 19.0, to which 
30 objects were classed as stars, and all fainter detections 
were assumed to be galaxies. We estimate the resolution 
of the ISAAC data as the mean Gaussian FWHM of the 
non-saturated stars, 0.46 arcsec. At this depth, the slope of 
the Galactic star counts is only ~ 0.1, (e.g. Roche et 
al. 1999) which by extrapolation would indicate the star- 
contamination at 19 < Ks < 22 to be ~ 63 stars, or 4.6 per 
cent of the galaxy sample. 

Figure 1 shows the differential number count of galaxies 
as a function of total Ks mag, shown with other counts 
and two models from Paper I. Saracco et al. (2001) also 
observed in the Ks band, and the 2.2/im results of Paper I 
and Moustakas et al. (1997) are plotted assuming a small 
colour term Ks — K ~ 0.04, estimated for a passive ERO 
model at 1 < < 2. 
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Figure 1. Differential number counts of galaxies against total 
Ks mag, for the ISAAC CDFS and ELAIS:N2 (Paper I) fields 
(shown well beyond completeness to illustrate the difference in 
depth), and from the even deeper surveys of MoustaJjas et al. 
(1997) and Saracco et al. (2001) (JC-band observations are plotted 
assuming Ks — K = 0.04). Also plotted are PLE (dotted) and 
R4, = Rm' = 0.3 merging (dashed) models (as in Paper I but 
with h = 0.7.) 



Our CDFS galaxy counts agree well with the Paper I 
ELAIS:N2 counts to Ks ~ 21, and reach ~ 1.5 mag deeper, 
turning over only at Kg > 22. Compaxison with the two even 
deeper surveys suggests we are near-complete to Ks ~ 21.5 
with some (~ 20 per cent) incompleteness at 21.5 < Ks < 
22. In most of the analyses of this paper we consider samples 
limited at total magnitudes of Ks = 22.0 and Kg = 21.5, for 
the more incompleteness-sensitive measures giving primary 
emphasis to the Ks < 21.5 sample. 



4 SELECTION OF THE ERO SAMPLE 
4.1 Identifying the EROs 

EROs are essentially galaxies with the red colours expected 
for high-redshift ellipticals. In Paper I they were selected on 
the basis of a colour R — K > 5.0, while some authors have 
used slightly stricter criteria of I — K > 4.0 or R — K > 5.3 
(e.g. Mannucci et al. 2002). Throughout this paper we select 
EROs as 1775—Ks > 3.92, which corresponds to our evolving 
E/SO model at 2 > 0.99 and is approximately equivalent (on 
the basis of a passive ERO model) to the R — K > 5.0 of 
Paper I and to I — Ks > 3.75 for the slightly longer 7-band 
(eg. as on the VLT FORSl instrument). 

The flrst step in identifying the EROs was to match 
each object in the Ks source list with its closest counter- 
part on the ACS 7775 image, using the RA and Dec co- 
ordinates given by the astrometry supplied with these im- 
ages. We allowed ofTsets up to a maximum 1.0 arcsec radius. 
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Matching the Ks source hst to their H and J counterparts 
was more straightforward, as the use of double-image mode 
meant there were no positional offsets. 

ERO selection was initially performed on a subset of the 
data and the candidate EROs examined by eye to check on 
the reliability of the procedure. The selection method was 
then modified where necessary. 

A list of EROs, numbering 155, was then selected as all 
galaxies detected in both Kg and 1775 , with total magnitude 
Ks < 22.0 and colour (from the 2.0 arcscc diameter aperture 
magnitudes) /tts — Ks > 3.92. Fixed-aperture colours were 
found to be more reliable for ERO selection than colours 
from total magnitudes, the reason being that in a significant 
fraction of galaxies SExtractor measured total magnitudes 
within very different sized apertures for Ks and 7775, giving 
rise to errors and selection biases. The mean offset between 
the Ks and 7775 centroids was only 0.225 arcsec. With this 
positional accuracy, and a background density of about 0.05 
7775 detections arcsec"^, the probability of chance super- 
position producing a false match (closer than the true 7775 
counterpart) is only ~ (7r0.225^ x 0.05), 0.75 per cent, for 
each galaxy. 

In addition, 53 galaxies with Ks < 22.0 had no detec- 
tion within 1.0 arcsec on the ACS data. Many of these will 
also be EROs, but not all - some, for example, are confused 
with neighbouring brighter objects in the shorter passband, 
causing their 7775 detection centroid to be shifted outside 
of the matching radius. To 'filter' those 53 objects, wo used 
IRAF 'qphot' to measure their aperture (2.0 arcsec diameter) 
magnitudes on the ACS images at the co-ordinates positions 
given by their ISAAC detections, and then only accepted the 
objects with (7775 — Ks)apcrture > 3.92, which numbered 43. 
Visual chocking indicated that this method appeared to ex- 
clude the confused objects. 

Our ERO sample then consisted of a total of 155 + 
43 = 198 galaxies with Ks < 22.0, of which 179 were 
brighter than the estimated completeness limit of Ks = 21.5. 
As none of the brighter {Ks < 19.0) objects classed as 
Galactic stars had 7775 — Kg any redder than 3.0, star- 
contamination of the ERO sample will be zero or very small 
(< 2). 

4.2 Colours of EROs 

ERO colours correspond to passive galaxies at approxi- 
mately 2 > 1, where the the observer-frame i^-band and the 
R or 7-band will bracket the large break at 4000A rest-frame. 
However, any type of galaxy, including a young starburst, 
could be this red if sufficiently dusty. 

Figure 2a shows models of I — Ks against redshift, which 
start forming stars 12.55 Gyr ago, at z = 6. We model galaxy 
spectral energy distributions using a version of the 'Pegase' 
package (Fioc and Rocca-Volmerange 1997; see Paper I). 
The elliptical galaxy model (solid line) forms all its stars in 
a 1 Gyr burst at 6 > 2: > 3.4, and the SO model with a 2 
Gyr burst continuing to z = 2.6. Both have E{B — V) = 

0. 4 mag dust extinction during the bursts and subsequent 
unreddened passive evolution, and show ERO colours at z > 

1. This choice of formation redshift appears reasonable in 
view of the discovery of star-forming galaxies in all intervals 
of redshift out to the current maximum of z = 5.78 (Bunker 
et al. 2003). 
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Figure 2. (a) Observer-frame 7775 — Ks colour (Vega system) 
against redshift for models representing evolving elliptical (solid) 
and SO (dashed lines) galaxies and dust-reddened young star- 
bursts (dotted/dot-dash lines), (b) Observed 7775 — Kg (in 2.0 
arcsec diameter apertures) for all Kg < 22 galaxies, with larger 
symbols for the 198 EROs. 

Dusty starburst models, represented as non-evolving 
starbursts observed at age 50 Myr, and heavily dust- 
reddened by £^(7? — V) = 0.8 and 1.0 mag, are also plotted, 
as dotted and dot-dash lines. These produce similar I — K 
colours to passive galaxies, especially at 2 > 1, and if very 
dusty they could have Irrs — Ks > 3.92 at 2 --^ 0.5-1.0. Older 
star-forming galaxies, or recent post-starbursts, would have 
intermediate spectra and could be EROs with more moder- 
ate dust-reddening. 

Figure 2b shows 7775 — Ks for all Ks < 22 galaxies, with 
larger symbols for the EROs, which, as in Paper I, appear at 
Ks ~ 18. Most 7775 < 25.5 EROs are detected on the ACS 
data, most at 7775 > 26 are not. 

4.3 Colour classification 

It is possible, with the H — Ks and J — Ks colours provided 

by the ISAAC data, to achieve an approximate separation 
of passive and dusty, star-forming EROs (pEROs and ds- 
fEROs). Pozzetti and Mannucci (2000) describe how these 
types can be separated on a plot of 7 — TsT against J — Ks, 
although this is only possible at 1 < 2 < 2, where the red- 
shifted 4OOOA break lies between the 7 and J bands. Here, 
we use essentially the same method, but also incorporate the 
77 — Ks colour, which extends this method to 2 ~ 2.5. 

Figure 3 shows 1775— Kg against J+H—2Ks (ie. the sum 
of J — Ks and H — Kg) with the diagonal line showing the 
adopted pERO/dsfERO divide, at J+H-2Kg = 0.86(7775- 
7^,) - 0.75. 

The majority of these EROs have colours consistent 
with either the passive models at 1 < 2 < 2.5 or the 



, ERO (I-deLected, I,,5-K,>3.92) 
3 ERO (not detected in I) 
„ other galaxy 
a star 




18 
Kg mag 
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,ERO (I„5-K,>3.92) 




, .. ..... *.* ■• •= " 



K<20 30<.-K<22 

. • • P^l^" 

' ■ •,- » . cisfl HO 

^Chandra source 



53,1 
R.A. (deg) 



Figure 3. Plot of Irrs-Ks vs. J+H-2Ks for all K < 22.0 EROs. 
Also plotted, with redshifts marked, are the elliptical (solid) and 
SO (short-dashed) models, with heavier lines indicating the red- 
shift range at which they are forming stars, and starburst models 
with reddening of E{B — V) = 0.8 (dotted) and 1.0 magnitudes 
(dot-dash) . The long-dashed diagonal line shows the adopted sep- 
arator of pERO and dsfERO types. 



E{B - V) ~ 0.8 starburst at z < 2. Using the plotted 
dividing line, we classify 109 EROs as pEROs and 89 as 
dsfEROs. It must be emphasised that these classifications 
are subject to large random error, and many EROs lie close 
to the dividing line, whether from scatter in the magnitudes 
or to genuinely intermediate spectra. To estimate the effec- 
tiveness of the classification we (using the same method as 
Mannucci et al. 2002) calculate for each ERO, the horizon- 
tal distance (i.e. A( J + H - 2Ks) from the pERO /dsfERO 
dividing line on Figure 3, and divide this by the statistical 
error in the EROs {J + H - 2Ks) colour. 

Many EROs are found to have strong, ~ 3-6a clas- 
sifications, but 47 classed as pEROs and 55 dsfEROs are 
within la of the lino, and hence about half of the individual 
classifications (102) must be regarded as weak. The result- 
ing uncertainties in the total number of dsfEROs would be 
smaller by ~ y'A'' and if we add this to the ^JN ol the total 
number, we estimate there to be 89±v^(55-|-89) = 89±12 ds- 
fEROs. There would of course be an additional uncertainty 
from the model-dependence of the dividing lino. However, 
our estimated dsfERO fraction, 45 ± 6 per cent, is consis- 
tent with previous estimates of 33-67 per cent from direct 
spectroscopy of Ks < 19.2 EROs (Cimatti et al. 2002a), 40- 
60 per cent from multicolour photometry of K' < 20 EROs 
(Mannucci et al. 2002), and 45-60 per cent from combined 
radio observations and colours of K' < 20.5 EROs (Smail et 
al. 2002). 



- « 



ft 



K<20 20<K<22 
9 » pERO 

i . dsfERO 

^Chandra source 



53.15 53.1 53.05 

R.A. (deg) 



Figure 4. (a) Distribution of colour-classifed pEROs and ds- 
fEROs, other detected galaxies and Chandra X-ray sources, on 
the area of the CDFS covered by the ISAAC data. The dashed 
circle indicates a possible cluster of EROs centered on the Chan- 
dra source XID:58 (see Section 7.1). (b) As (a) but with the field 
boundaries shown and the non-ERO galaxies omitted, for clarity. 



4.4 Distribution on the Sky- 
Figure 4 shows the distribution of the ERO sample, and the 
73 Chandra detections (Giacconi et al. 2002) on the ISAAC 
mosaic area of the CDFS. The EROs appear somewhat clus- 
tered on small scales, and a particularly obvious overdensity, 
centred on Chandra source XID:58 (R.A. 3''32'"11.85'' Dec - 
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Table 1. Observed number counts (number of galaxies Ng and 
surface density p (in units deg~^mag~^ , with ^JN errors) of EROs 
(7775 -Ks> 3.92 galaxies) on the CDFS. 







p 


< 17.0 








17.0-17.5 


1 


143 ± 143 


17.5-18.0 


1 


143 ± 143 


18.0-18.5 


4 


572 ± 286 


18.5-19.0 


18 


2573 ± 606 


19.0-19.5 


20 


2859 ± 639 


19.5-20.0 


25 


3574 ± 715 


20.0-20.5 


28 


4003 ± 756 


20.5-21.0 


40 


5718 ± 904 


21.0-21.5 


41 


5861 ± 915 



27:46:29.14), is highlighted on the plot and discussed further 
in Section 7.1. 



5 NUMBER COUNTS OF EROS 

5.1 Observed counts on the CDFS 

The differential number counts for EROs on the CDFS are 
given in Tabic 1 and plotted on Figure 5. Firstly, the three 
ERO counts plotted here are consistent within the error 
bars. Our CDFS counts are closer to the uncorrected count 
of Paper I, suggesting the incompleteness corrections may 
have been overestimated. For this paper, rather than at- 
tempting to correct for incompleteness, we have assumed 
completeness to Kg = 21.5, and disregarded the final point 
in fitting slopes and models (incompleteness will appear 
more suddenly for this data due to its more uniform depth) . 

The most obvious feature on the plot is a marked flat- 
tening of the count slope, seen in all three sets of ob- 
servations but most apparent with the CDFS. The count 
slope best-fitting the the ELAIS:N2 and CDFS counts at 
17.0 < Ks < 19.5 is 7 = 0.59 ± 0.11, and flattens to 
7 = 0.16 ± 0.05 at 19.0 < < 21.5 (fit to CDFS counts 
only). 



5.2 ERO counts: Models and interpretation 

In Paper I we presented simple models of the number counts 
of EROs, based on the assumption that both pEROs and ds- 
fEROs are progenitors of the present-day E/SO galaxies, and 
hence that their luminosity function (LF) at any redshift can 
be linked to the LF of local E/SOs by means of a particular 
evolutionary model. We adopted the local E/SO LF derived 
by Kochanek et al. (2001) in the K-h&nA, from the 2MASS 
survey, a Schechter function with = —24.31, a = —0.92 
and 0* = 0.0015435 Mpc~3 for h = 0.7. 

The ERO counts were modelled by evolving this LF 
with the galaxies split 50:50 between our elliptical and SO 
models of passive L* evolution (see Section 4.2), and includ- 
ing only the galaxies within the redshift ranges in which their 
respective models give ERO colours (i.e. approx. z > 1). 
In the Pure Luminosity Evolution (PLE) model only L* is 
evolved, while </>* and a remain constant from the formation 
redshift onwards. The post-starburst fading of dsfEROs is 



1 I I , , , I I I 

ERO counts 

oRoche et al. (2002) (ELAIS:N2) (R-K>5.0) 

ARociie et a!. (2002) 'corrected' 

□ Daddi el al. (2000) (R-K,>5.0) 

.This work (CDFS) (1„5-K.>3.92) . ' ' 



'7 
L 



Figure 5. A'-band differential number counts for EROs selected 
as /776 — Kg > 3.92 on the CDFS and at comparable limits of 
R - > 5.0 on the ELAIS:N2 field (Paper I/Roche et al. 2002), 
and R—Ks > 5.0 by Daddi et al. (2000). These axe compared with 
the counts of Ij-rs—Ks > 3.92 galaxies predicted by PLE (dotted), 
non-evolving (long-dashed), merging with = Rm = 0.3 (short- 
dashed), and 'M-DE' merging with Rm = 0.3 R^ = -0.49 (dot- 
dashed) models. Error bars are ^N. 



assumed to be similar in degree to the passive evolution of 
the pEROs. 

In Paper I, the ELAIS:N2 ERO counts were significantly 
lower than expected if all present-day E/SOs evolved by this 
PLE model, and much closer to a non-evolving model. How- 
ever, the compact angular sizes of most EROs implied sig- 
nificant evolution in their surface brightness. 

A merging model was also considered, which combined 
the PLE model's evolution of luminosity (per unit mass) 
with an increase with redshift in number density (</>*) and 
a corresponding decrease in mass (parameterized as J?^ = 
Rm = 0.3; see Paper I) of the galaxies. The merger rate 
and its evolution were based on the observational estimate 
of Patton et al. (2002). This model also overpredicted the 
ERO counts. However, by keeping the merger model's mass 
evolution at Rm ~ 0.3 and varying the density evolution, 
a reasonable best-fit to the ERO counts on ELAIS:N2 was 
obtained for R^ = -0.46(±0.10). 

In this model, described as 'merging and negative den- 
sity evolution' (M-DE), the comoving number density of pas- 
sive galaxies (or any galaxies of similar colour) gradually 
decreases with redshift. Physically, this could result from 
some fraction of the present-day massive ellipticals forming 
at high (e.g. z > 3) redshifts, with the remainder forming 
from bluer galaxies (e.g. through interactions and mergers) 
at all intermediate redshifts. 

With the deeper data available for the CDFS, we should 
now be able to improve the constraints on ERO evolution. 
We recompute the models of Paper I for a ft. = 0.7 Universe 
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and ERO selection as hrs - Ks > 3.92. The M-DE model 
best-fitting the CDFS counts has a density evolution param- 
eter of = — 0.49(±0.06) (almost identical to the best-fit 
model of Paper I). 

Figure 5 compares the new ERO counts with PLE, 
merging {R^ = Rm = 0.3), non-evolving and M-DE {R^ = 
—0.49) models. As in Paper I, the PLE and merging mod- 
els overpredict the observations. The deep CDFS counts arc 
consistent with the non-evolving model at brighter magni- 
tudes, but fall below it at Kg > 20. This suggests that there 
is a real decrease in the comoving 0* of ERO gala^xies with 
redshift, and not only the absence of strong luminosity evo- 
lution. 

Of these models, M-DE is the only one consistent with 
the plotted counts. In the M-DE model the evolution of the 
characteristic mass (from merging) is (as in Paper I), m*(z = 
1,2,3) = (0.714, 0.558, 0.493)m*(« = 0), and the comoving 
number density of passive/rod galaxies, with R^ = —0.49, 
evolves as ^'(z = 1,2,3) = (0.577, 0.386, 0.315)m* (a = 0). 



6 CLUSTERING OF THE EROS 
6.1 Calculating u){9) 

To investigate the clustering of EROs on the CDFS, and 
compare it with other galaxies, we calculate the angular cor- 
relation function, cij{6), for (i) all galaxies (of any colour), 
(ii) all EROs, and (iii) pEROs and dsfEROs considered sepa- 
rately. The method adopted is exactly as described in Paper 
I, with the use of Nr = 50000 random points. 

In addition, (iv) we calculate the angular cross- 
correlation of the pEROs and dsfEROs, as 



where Npd(Oi) is the total number of pERO-dsfERO pairs 
(centering on the pEROs) in the 9i interval of angular sep- 
aration, Npd(9i) the number of pERO-random point pairs, 
Nd the number of dsfEROs and Nr the number of random 
points. As an additional check, we also evaluated the cross- 
correlations with centering on the dsfEROs, i.e. as 



<jJdp{di 



NdpiOi) Np 

NdriOi) Nr 



- 1 



which gave essentially the same results - LOdp — oJpd at all 0, 
to within a small fraction of the statistical errors, ~ O.lcr. 

The 3-dimensional two-point correlation function of 
galaxies is an approximate power-law in physical separa- 
tion (r), C(r-) = (7^)"^ ® (Peebles 1980), at separations of 
up to a few Mpc. following from this, the U3{d) of galax- 
ies will also be an approximate power-law uj{9) ~ Au9~°'^, 
at the relatively small 6 considered here (< 0.1 dcg). We 
therefore (as in Paper I) express the observed io{9) as an 
amplitude A^. at 6 — 1 dcg, obtained by fitting with the 
function 'A„(6'"°-* - 12.24)' (where 12.24 is the 'integral 
constraint' as determined for the ISAAC field area; see Pa- 
per I), over the range of separations 2 < 9 < 32 arcscc. Error 
bars were estimated by the method described in Paper I, us- 
ing 12 'sub-areas', and are generally a little larger than -^/N 
errors. 



I • • • • I ' 

.,K<ai,5 EROs 

iiK<22.0 EROs 

,K<21.5 all galaxies 



(b) 



^K<21.5 pEROs 
...K<21.5 dsfEROs 
-A cross— correlation 



1 1.5 

Separation, log d arcsec 



Figure 6. (a) Angular correlation functions i^{6) for EROs and 
galcixies of all colours on the ISAAC CDFS data, at the ap- 
proximate completeness limit of Kg = 21.5, and for EROs at 
the fainter Kg = 22.0 limit, with the dashed, dotted and dot- 
dash lines showing respective best-fitting functions of the form 
'Au,{0-°-^ - 12.24)'; (b) iv{0) of K < 21.5 EROs of pERO and 
dsfERO colour-classification, and their cross-correlation function, 
again with the dashed, dotted and dot-dash lines showing best- 
fitting 'A^{0-°-^ - 12.24)'. 



6.2 oj{9) Results 

Table 2 gives our A„ amplitudes for a range of magnitude 
limits (these results may be less reliable at Ks > 21.5 due to 
incompleteness), and Figure 6 shows u}{9) with fitted power- 
laws. Wo obtain statistically significant, ~ 3a detections of 
clustering, both for all gala^xies and for the EROs to Kg = 
21-22. 

Table 2 also gives the A^ for the pERO and dsfERO 
subsamples. Clustering is detected for pEROs and consis- 
tent with a similar amplitude to the full ERO sample, but 
is at best 2.3a. The detection of dsfERO clustering and the 
cross-correlation are even weaker, and would be consistent 
with this same amplitude or one only half as groat. Clearly, 
more data are needed before we can make a meaningful com- 
parison between the clustering properties of pEROs and ds- 
fEROs. 



6.3 Interpretation of the u>{9) 

The observed ui{9) of any sample of galaxies will depend 
on their intrinsic clustering in 3D space, described by the 
two-point correlation function ^(r), and their redshift dis- 
tribution N{z). If ^(r) is represented by the simple model 



i{r,z) = {r/ro)-\l + z) 



-(3-He) 



where ro is the local correlation radius, 7 ~ 1.8 (observation- 
ally) and e represents the clustering evolution (e = is stable 
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Figure 7. The scaling of hj(6') amplitude with K magnitude limit, 
for full i^-limited samples of galaxies from our ISAAC CDFS 
data, with previously published results, and models (all with 
ro = 5.85 h^QQ Mpc) of stable (e = 0) and comoving (e = —1-2) 
clustering and an intermediate (e = —0.4) evolution. 

and e = —1.2 is comoving clustering), then Limber's formula 
(see eg. Efstathiou et al. 1991) gives u{e) = ^u;6'"*^"^\ 
where 

Jo X' ^{z) Jo 

where x{z) is the proper distance and = 3.679 for 7 = 
1.8. 

Before considering the EROs we briefly investigate the 
io{9) scaling of full i^-limited galaxy samples. Figure 7 shows 
our observations for the CDFS, and previously published re- 
sults, which agree well. No correction has been applied to 
the CDFS amplitudes for star contamination, but this will 
be very small, e.g. A(log Au,) = 0.04 for the estimated num- 
ber of stars (Section 3.4). Also plotted are models based on 
the n{z) from our = Rm = 0.3 merging model (see Paper 
I), a local correlation radius of ro = 5.85 Mpc (from the 
7 < 22 survey of Cabanac, de Lapparent and Hickson 2000), 
and showing the —1.2 < e < range of clustering evolution. 
The observations appear most consistent with e ~ —0.4-0, 
with comoving clustering excluded. 

Figure 9 shows the Au, scaling of EROs selected as hrs — 
Ks > 3.92 or R- K > 5.0. Our results for the CDFS are 
consistent with the scaling trend seen in previous surveys. 

Any estimation of the correlation radius (ro) from io{9) 
will inevitably be dependent on the assumed model for 
N{z). We represent the ERO N{z) using the 'M-DE' model 
(Figure 8) Currently there is insufficient spectroscopic data 
for EROs to fully test the accuracy of this N(z). We can, 
howver, compare our model with the Cimatti et al. (2002b) 
redshift (mostly spectroscopic, some photometric) survey for 
all types of galaxy at an effective limit ~ 19.75. Figure 8 



Figure 8. N{z) for 7775 — Ks > 3.92 EROs as given by our 
= -0.49 'M-DE' model, at limits = 19.0, 19.75, 21.0 and 
22.0. Also shown for comparison (dotted histogram) is the z > 1 
tail of the Cimatti et al. (2002b) N{z) for all types of galaxy to 
an effective limit K ~ 19.75, normalized to the ERO model at 
this magnitude. 

shows the z > 1 tail (138 galaxies) from this N{z), which 
is not greatly different from our model and the two may 
be consistent if the fraction of EROs in the survey sam- 
ple is highest at z > 1.2. Our model N{z) is also in good 
agreement with the physically-based Granato et al. (2003) 
prediction for spheroidal galaxies only. 

Ifence, assuming our model to be a good representation 
of the ERO N{z), we vary the ro normalization to best- 
fit (minimize x^) all the plotted data points. We estimate 
ro for EROs to be 12.5(±1.2) h'^ Mpc (^ = 1.04) for 

comoving and ro = 21.4(±2.0) Mpc (^ = 0.77) for 
stable clustering. 



7 ERO CLUSTERS AND PAIRS 
7.1 A possible cluster of EROs 

In Paper I we noted a possible cluster of EROs in the 
ELAIS:N2 field, centered on a bright ERO which is also 
a Chandra source and radio galaxy. 

In the CDFS, we find another overdensity of EROs, cen- 
tered on Chandra source XID:58 (R.A. 3''32'"11.85'' Dec - 
27:46:29.14), a somewhat fainter ERO {K^ = 20.94, H = 
21.66, J = 22.92 and 1775 = 24.97). Within a 20 arcsec 
radius of XID:58 there are a total of 10 K < 22.0 EROs 
compared with 1.37 expected from their mean surface den- 
sity on the field. There is no accompanying overdensity of 
bluer galaxies, with 9 seen compared to 9.83 expected. 
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Table 2. Galaxy ui{6) amplitudes A,^ (in units of 10~^ at one degree) for full Ks-limited samples of galaxies, EROs (7775 — Kg > 3.92), 
EROs classed as pEROs and dsfEROs, and the pERO/dsfERO cross-correlation, to a series of magnitude limits. Ng is the number of 
galaxies in each sample. The A^ of EROs is not given for the first two limits as the samples are too small for any detection of clustering. 





All galaxies 




All EROs 




pEROs 




dsfEROs 


pERO X dsfERO 


limit 






Ng 




Ng 




Ng 


A.QJ 


A.QJ 


20.00 


467 


12.53 d= 4.89 


69 




41 


28 








20.50 


638 


10.56 ± 3.23 


97 




56 


41 








21.00 


898 


8.92 ± 2.28 


137 


53.5 ± 19.7 


79 


79.4 ±48.7 


58 


72.9 ±34.9 


31.3 ±30.0 


21.50 


1233 


6.98 ± 2.50 


179 


60.3 ± 19.0 


101 


56.0 ± 24.2 


78 


43.0 ±31.0 


39.6 ±21.2 


22.00 


1619 


5.05 ±2.01 


198 


48.9 ± 15.8 


109 


45.0 ± 19.6 


89 


23.5 ± 22.5 


31.2 ± 23.4 



Table 3. EROs detected as X-ray sources, with X-ray properties from Giacconi ct al. (200) (: Chandra source XID number, RA and 
Dec of the corresponding ERO on the Kg image. Kg magnitude (total) , I775 — Kg colour (aperture) , pERO /dsfERO colour classification 
(brackets denote weak, < la classifications). Soft and Hard X-ray fluxes in units of 10~^® erg s~^cm~^ (or \a upper limits), and X-ray 
hardness ratio HR (see text). 



XID 




R.A 




Dec. 






Ks 




,^77,5 - 


Ks 


p/dsf 


F(0. 5-2.0 keV) 


F(2-10 keV) 


HR 


26 


3: 


:32: 


:39, 


.75 


-27:46: 


:11, 


.25 


19, 


.64 ±0, 


,09 


4.05 ± 


0, 


,11 


(P) 


7.35 ±0.750 


28.9 ±4.07 


-0.23 ±0.08 


58 


3: 


:32: 


:11, 


.76 


-27:46: 


:28, 


.14 


20 


.94 ±0, 


17 


4.04 ± 


0, 


,20 


(dsf) 


6.36 ±0.668 


18.8 ±3.26 


-0.36 ±0.09 


79 


3: 


:32: 


:38, 


.04 


-27:46: 


:26, 


.26 


20, 


.92 ± 0, 


.11 


4.83 ± 


0, 


,22 


(P) 


8.06 ± 0.785 


20.7 ±3.64 


-0.42 ±0.08 


86 


3: 


;32: 


;33, 


.84 


-27:45: 


;20, 


.45 


21, 


.22 ±0. 


.16 


4.82 ± 


0, 


.34 


(P) 


1.22 ±0.345 


6.96 ± 2.43 


-0.05 ±0.22 


108 


3: 


;32: 


;05, 


.76 


-27:44: 


;46, 


.91 


21, 


.28 ± 0. 


.12 


4.16 ± 


0, 


.20 


P 


2.94 ± 0.513 


8.73 ± 2.92 


-0.36 ±0.17 


153 


3: 


;32: 


;18, 


.32 


-27:50: 


;55, 


.25 


18 


.58 ± 0, 


,04 


4.34 ± 


0, 


.05 


P 


1.88 ± 0.403 


46.6 ± 4.43 


0.59 ± 0.07 


188 


3: 


;32: 


:22, 


.53 


-24:49: 


:49, 


.75 


18 


.17 ±0, 


,04 


3.94 ± 


0, 


,04 


dsf 


< 0.616 


6.43 ±2.37 


1.00-0.81 


221 


3: 


:32: 


:08, 


.91 


-27:44: 


:24, 


.81 


21 


.76 ± 0, 


,16 


5.63 ± 


0, 


,84 


(dsf) 


1.09 ±0.360 


< 5.12 


-1.00+0-84 


253 


3: 


:32: 


:20, 


.05 


-27:44: 


:47, 


.20 


20 


.00 ± 0, 


,09 


3.97 ± 


0, 


,11 


(dsf) 


0.829 ± 0.320 


26.0 ±3.52 


0.66 ±0.11 


254 


3: 


:32: 


:19, 


.80 


-27:45: 


:18, 


.72 


21 


.48 ± 0, 


,21 


4.06 ± 


0, 


,28 


(dsf) 


< 0.573 


18.4 ±3.06 


1.00-0.32 


513 


3: 


:32: 


:34, 


.01 


-27:48: 


:59, 


.74 


21 


.19 ± 0, 


,17 


4.49 ± 


0, 


,25 


dsf 


0.847 ±0.331 


8.60 ±2.60 


0.23 ±0.23 


515 


3: 


;32: 


;32, 


.17 


-27:46: 


:51, 


.49 


21, 


.99 ± 0, 


.19 


4.10 ± 


0, 


.33 


(dsf) 


0.843 ± 0.304 


12.7 ±2.74 


0.41 ± 0.17 


572 


3: 


;32: 


;22, 


.13 


-27:48: 


;11, 


.47 


21 


.17 ±0. 


,16 


4.92 ± 


0, 


.29 


(dsf) 


1.00 ± 0.368 


< 5.39 


-1.00+0-90 


574 


3: 


;32: 


;31, 


.53 


-27:48: 


;53, 


.81 


20 


.25 ±0. 


.11 


4.28 ± 


0, 


.13 


(dsf) 


0.787 ±0.319 


< 4.44 


_100+0-93 


600 


3: 


:32: 


:13 


.81 


-27:45: 


:25, 


.63 


18, 


.66 ± 0, 


,03 


4.20 ± 


0, 


,04 


P 


< 0.548 


18.8 ±3.26 


1.00-0.66 


609 


3: 


:,32: 


:36 


.17 


-27:50: 


:,36 


.93 


19 


.64 ± 0, 


,07 


4.49 ± 


0, 


,09 


(P) 


< 0.687 


9.31 ± 2.63 


1.00_„.65 


623 


3: 


:32: 


:28, 


.51 


-27:46: 


:58, 


.17 


20 


.60 ± 0, 


,09 


4.51 ± 


0, 


,17 


dsf 


0.809 ± 0.300 


< 4.49 


-1.00+0-92 
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Figure 9. The scaling o{ ll>(9) amplitude with K magnitude limit 
for EROs selected as hrsl - > 3.92 or R-K > 5.0, from our 
ISAAC CDFS data and from previously published results. These 
are compared with M-DE models with either stable (e = 0) or 
comoving (e = —1.2) clustering and ro best-fitted to the plotted 
data points. 
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Figure 10 shows this region on the ISAAC image with 
EROs and Chandra sources marked, and Figure 11 shows 
these 10 EROs on a colour-colour plot. 

A recent photometric redshift estimate for the XID:58 
galaxy (Mainicri ct al., PhD thesis, in prop) gives z ~ 1.44. 
If the ovcrdensity is at this redshift, a 20 arcsec radius is 
0.233 Mpc for h = 0.7. On the colour-colour plot, six of the 
EROs within this radius appear clustered around a passive 
model at a ~ 1.5. It seems very likely that these pEROs and 
the X-ray source form a single cluster. 

Of the 3 dsfEROs, one, the = 20.39 galaxy, has 
colours at least consistent with membership of this cluster. 
The colours of the other two suggest lower z, foreground 
galaxies. A second X-ray source within the cluster area, 
XID:149, has a spectroscopic redshift of only 1.033 (Szokoly 
et al. 2003). 

On Figure 11 the XID:58 galaxy is classed as a ds- 

fERO but effectively lies on (O.OIct from) the adopted 
pERO/dsfERO divide. A dsfERO classification seems more 
likely on the basis of (i) a ~ 0.5 mag excess in the /775-band 
{\rest — 320nm) relative to the E/SO model colours at its 
'photo-z', z = 1.44, and (ii) the apparent morphology on 
the ACS image which docs not resemble a giant elliptical, 
but rather is asymmetric and very compact. It z = 1.44 is 
assumed to be its correct redshift, we estimate the i?-band 
absolute magnitude as Mr ~ —22.2 (for h = 0.7), from 
the observed magnitude closest to restframe R, H = 21.66, 
and our passive model. This would, in our M-DE model, be 
consistent with an L ^ L* galaxy at this redshift. 

The X-ray flux, F(0.5-2.0 keV) = 6.36(±0.668) x 10""' 
and F(2.0-10.0 keV) = 1.88(±0.326) x 10"^^ ergs cm'^s"^ 



Figure 10. Greyscale plot of the 40 x 40 arcmin region of the 
ISAAC Ks image centred on the Chandra source XID;58, on 
which X-ray source centroids are labelled (diagonal crosses) with 
XID numbers from Giacconi et al. (2002) and EROs (circles) by 
pERO/dsfERO classification and Kg magnitude (in brackets). 



would, for a /j, oc v ^ SED, corresponds to Lx(0.5- 
10.0 koV) ~ lO''-'^-''' ergs s"^ for h = 0.7, which is too high 
for a non-AGN starburst or elliptical (see Section 8). The 
X-ray emission cannot be from hot intracluster gas as it is 
consistent with a point-source in the Chandra data (Giac- 
coni et al. 2002). The non-detection of a diffuse component 
might not be unexpected at z ^ 1.5, as the X-ray luminosity 
function of clusters is known to undergo rapid negative evo- 
lution with increasing redshift (e.g. Henry 2003). The high 
X-ray and modest optical luminosity of the central ERO, 
combined with the morphology and colours all suggest it 
is a dusty starburst or post-starburst galaxy (likely a recent 
post-merger) , and the host of an obscured AGN - the source 
of the X-ray emission. 

In summary, wo find evidence for the existence of a 
gala^xy cluster at 2; ~ 1.5, containing 7 or 8 EROs with 
K < 22, including an X-ray luminous obscured AGN (Chan- 
dra source XID:58) and several passive galaxies. A notable 
difference from the possible z ~ 1.1 cluster reported in Pa- 
per I is that the X-ray source does not lie in the central giant 
elliptical but a in a relatively small, irregular and probably 
dsfERO type galaxy. Some of the cluster ellipticals are up 
to 2 mag brighter than XID:58 and would be within the 
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Figure 11. Colour-colour plot of the 10 EROs within a 20 arc- 
sec radius of XID;58, each labelled by Ks magnitude (as on 
Figure 10). The long-dashed diagonal line shows the adopted 
pERO/dsfERO divide, and the horizontal long-dashed line the 
colour limit for ERO selection. Redshift loci for the E and SO 
models are plotted as solid and short-dash lines, with the bolder 
lines indicating their high-z starburst phases. Redshift loci for the 
dusty starburst models are plotted as dotted (for E(B — V) = 0.8) 
and dot-dash (for E{B — V) = 1.0) lines. Redshifts are marked 
on each model. 



reach of spectroscopy, allowing the existence of a cluster to 
be verified. 



7.2 Close pairs of Galaxies, and the Merger rate 

The number of pairs of galaxies at very close separations, 
< 5 arcsec (~ 40 h^^ kpc at 1 < a < 3) is related primar- 
ily to the frequency of interactions and mergers. Patton et 
al. (1992) estimated that for an optically selected sample 
of galaxies at « ~ 0.3, there were 0.0321 ± 0.007 interact- 
ing companions per galaxy, and that this fraction and the 
merger rate were evolving (at 2; < 1) as oc {l + zf-''^. This 
result is incorporated in our merging and M-DE model as a 
merger rate Rm.{z) = 0.3(1 -|- z)^'^ per Hubble time. It also 
provides a useful point of comparison for the EROs. 

We investigate the close-pair counts within the ERO 
sample, using a method originally described by Woods, 
Fahlman and Richer (1995). For each ERO-ERO pair, the 
probability P of occurring by chance (in a random distribu- 
tion) is estimated as 



P = 1 — exp[— p(< m2)fi,3e] 



(1) 



where p(< m2) is the surface density of galaxies in the sam- 
ple brighter than m2, the apparent magnitude of the fainter 
galaxy of the pair, 

Qpe is the area of the annulus around the brighter 
galaxy between /?, an angular separation cut-off below which 



Figure 12. Histogram of the number of close ERO-ERO pairs 
with P < 0.05 (as defined in the text) as a function of pair sep- 
aration, as calculated for our 179 EROs with Kg < 21.5 and 
/ — Kg > 3.92. The plotted points show the simulation prediction 
for the same number of galaxies distributed randomly (where the 
arror bars show the error for a single sample of 179) . 



individual objects cannot be resolved (/3 ~ 1 arcsec here), 
and 9, the pair separation ( Qpg = tv6^ — 1x0^ away from any 
edges or holes). We estimate Q-^e numerically by counting 
points distributed randomly over the observed field area. 

Galaxy pairs with P < 0.05 arc then counted as close 
pairs with a high probability that they are interacting. As 
faint companions of brighter galaxies are especially likely to 
be missed near the detection threshold, this analysis is lim- 
ited to the Ks < 21.5 EROs. Amongst these 179, we classify 
20 close pairs, of which 9 are pERO-pERO, 3 are dsfERO- 
dsfERO and 8 are mixed-type - there is no indication that 
dsfEROs are more likely than pEROs to be in interacting 
pairs. 

To estimate the number of close pairs resulting from 
chance coincidence, the same analysis is performed on 100 
simulations with the same number of galaxies and magni- 
tude distribution as the real sample, but randomized posi- 
tions. Averaging over all simulations, we find a random dis- 
tribution has 8.98±3.36 pairs with P < 0.05 (where the error 
bar is for a single simulation). Figure 12 shows the observed 
and random-simulation P < 0.05 pair counts for Kg < 21.5 
EROs - the observed counts show an excess at 1-5 arcsec. 
Subtracting the random expectation from the observed pair 
count gives 11.02 ± 3.36, and hence the rmmber of interact- 

-txae ^ o.062±0.019. 



ing ERO companions per ERO as 
This is higher than the Patton et al. (2002) estimate for 
optically-selected galaxies at 2: ~ 0.3, but with (1 -I- z)'^'^ 
evolution, this would reach 0.062 at z — 0.73. 

However, the P < 0.05 pair count is probably under- 
estimated due to image merging causing some of the very 
closest {0 < 2 arcsec) pairs to be missed. For a rough esti- 
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mate of this effect, we can assume the true ratio of galaxy 
pairs to the random prediction at ^ < 2 arcsec is no lower 
than at 2 < < 4 arcsec, where observationally it is 
At < 2 arcsec, we counted 2 pairs compared to the ran- 
dom prediction of 1.08, but from the above ratio wc should 



have 



X 1.08 = 3.96, an additional 1.96, increasing the 



paired fraction to = 0.073 ± 0.019 compan- 

ions per ERO. The Patton et al. (2002) evolution reaches 
this figure at « = 0.86, which is still less than the typical 
ERO redshifts. 

On this bEisis we tentatively conclude that the incidence 
of visibly interacting or merging pairs amongst EROs is com- 
parable to, and not significantly greater than, that occurring 
in an optically-selected sample of galaxies at the same z >1 
redshifts. 



8 CORRELATION OF EROS WITH X-RAY 
SOURCES 

The CDFS was surveyed with Chandra in two passbands, 
0.5-2.0 keV ('Soft') and 2-10 keV ('Hard'), reaching (in a 
total of 942 ksoc exposure time) respective source-detection 
limits (for a photon index F = 1.375, the mean for the de- 
tected sources) of 5.5 X 10"^^ and 4.5 x 10"^^ erg s ^cm 
(which are 12 and 15 counts Msec~^). Giacconi et al. (2002) 
catalog the X-ray sources, numbering a total of 346 detec- 
tions in the Soft (307) and/or Hard (251) bands, of which 
73 are in the area of our ISAAC data. For each source, an 
X-ray hardness ratio HR is given as where H and S 
are the photon counts in the Hard and Soft bands, and can 
take any value from HR = — 1 (detected only in Soft band) 
to HR = +1 (detected only in Hard). 



8.1 Calculating the cross-correlation 

We cross-correlate these 73 X-ray sources with the EROs 
to (i) identify which EROs are also X-ray sources and (ii) 
investigate whether there is any association between EROs 
and X-ray sources at non-zero separations. 

The high spatial resolution of Chandra, FWHM ~ 0.7 
arcsec near the aim point, provides source positions suffi- 
ciently accurate to be matched to individual galaxies, even 
to Ks = 22. When the Chandra sources wore correlated with 
the i^s-band source list, we found || to have probable Kg- 
band IDs within < 2.5 arcsec. We also found a systematic 
offset of ~ 1.3 arcsec in the Giacconi et al. (2002) X-ray cen- 
troids. To remove this, and optimize the positional match- 
ing, wo fitted an astrometric transform between these X-ray 
centroids (for 65 sources) and the positions in pixels of their 
closest Ks-h&nd counterparts. The rms scatter in this trans- 
form was only 0.44 arcsec. 

The Chandra source positions could then be trans- 
formed to accurate pixel co-ordinates on the Ks image. The 
source-ERO cross-correlation function (see Section 6.1) was 
evaluated as 

where Nxe{Oi) is the total number of Chandra-ERO pairs 
(centering on the Chandra sources) in the 6i interval of an- 
gular separation, Nxr{9i) the number of C/iandra-random 



pairs, Ne the number of EROs and Nr the number of ran- 
dom points (50000 here). 

8.2 X-ray detected EROs 

EROs were identified directly with Chandra sources where 
the Ks band and X-ray positions coincided within 2.0 arcsec. 

A total of 17 direct matches were found (each of these was 
checked by eye to verify selection of the correct counterpart), 
of which only 1.0 would (on the basis of the mean surface 
density of EROs) be expected from chance coincidence. 

If there are 17 matches, and ~ y| = 94 per cent are 
genuine, this means that 8.1 ± 2.1 per cent of the Ks < 
22.0 EROs are X-ray sources above the Chandra 942 ksec 
detection limit. Of the 17, we classed 7 as pEROs and 10 as 
dsfERO, and so estimate the X-ray detected fractions of each 
type as 6.4 ± 2.4 and 11.2 ± 3.6 per cent respectively. This 
is higher for dsfEROs but the difference is not statistically 
significant. These fractions are not significantly different at 
Ks < 21.5 (to which there are 14 matched EROs). Our result 
also means that — 



— = 23.5 ± 5.9 per cent of the Chandra 
sources detectable on the Ks image exhibit ERO colours, 
whether from a passive or a dusty host galaxy. 

Table 3 gives the X-ray and optical ID properties of 
these 17 EROs, and Figure 13 shows total 0.5-10.0 keV 
fiuxes (i.e. the sum of the 0.5-2.0 and 2.0-10.0 keV fluxes, 
assuming a zero fiux for any band in which there is no detec- 
tion) against HR. This plot can be compared directly with 
Figure 4 of Alexander et al. (2002) and Figure 2 of VignaU 
et al. (2003). 

The total X-ray flux and HR can help to classify the 
source as an AGN, starburst or elliptical galaxy. A flat spec- 
trum (photon index F — 1.0) is HR ~ —0.11. Both unob- 
scured AGN and starburst galaxies typically have ocv'^ 
(F = 2.0) SEDs, which is HR ~ -0.35. Obscured AGN will 
have a harder HR, depending on the absorption column, 
Nh and redshift. Starbursts have lower X-ray luminosities 
Lx < 10*^ ergs s~^, and hence at « > 1, fluxes F(0.5- 
10 keV) < 4 X 10"^® erg s"^cm"^, while AGN tend to be 
more luminous. Passive ellipticals, containing hot gas, have 
a similarly soft HR to the starbursts and even lower Lx ■ 

Of these 17 ERO X-ray sources, 13 (aU 7 pEROs, and 
6 dsfEROs) are detected in the Hard X-ray band (9 in both 
bands and 4 in Hard only). These 13 have total fluxes _F(0.5- 
10.0 keV) > 6 X 10"^'' ergs s^^cm"^, sufficient to be in- 
dicative of AGN, which must be obscured to some degree 
for the galaxies to exhibit ERO colours. Their HRs show 
a similar wide spread to the X-ray detected ERO sample 
of Alexander et al. (2002), and - as the fiux limits and 
hence (presumably) redshift ranges of these samples are sim- 
ilar correspond to the same range of absorbing columns, 
Nh ^ lO^'-'^-^^^-^cm-^ The pERO/dsfERO mixture sug- 
gests considerable diversity in the star formation histories 
of the AGN hosts. 

Alexander et al. (2003) have recently found that many 
(> 36 per cent of) bright Fgso^im > 5 mjy sub-mm sources 
are themselves detectable in the deepest (~ 1 MSec) X-ray 
surveys. The X-ray luminosities of these galaxies, Lx ~ lO''^ 
ergs s~^, indicate that at least a majority contain AGN, but 
are lower than QSO luminosities, and the sources also have 
high absorbing columns, Nh > 10^^ cm^^. These properties 
would be consistent with forming ellipticals, still in the pre- 
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Figure 13. Detected H + S X-ray flux (0.5-10.0 keV), against 
hardness ratio HR), for Chandra sources on the CDFS, with the 
17 identified with EROs indicated by large symbols and labelled 
with XID number. 

ejection, dusty phase. It seems possible that the 5 dsfEROs 
in our sample with HR > include similar proto-galaxies 
(the surface density of sub-mm sources detected with 1 MSec 
on Chandra corresponds to 5 on our sample area), but we 
require sub-mm data to verify this. 

We can, however, compare the X-ray properties of our 
EROs with the results of Alexander et al. (2002), who in 1 
Msec Chandra survey of the CDFN detected 6 of a total of 29 
HK' < 20 A 29 EROs {I - K > 4.0) in a 70.3 arcmin^ area. 
This HK' — 20.4 limit is approximately Ks = 20.1, and 5 of 
our X-ray detected EROs are brighter than this (Table 3). 
Hence, at an equivalent combination of X-ray and Ks flux 
limits, Alexander et al. (2002) and ourselves find consistent 
surface densities of EROs, = 0.085 ± 0.035 arcmin"^ 
and = 0.099 ± 0.044 arcmin^^ respectively. Alexander 
et al. (2002) also detected 1 of the 9 EROs in a deeper 
sample covering 5.3 arcmin^, and this is also consistent with 
our surface density at K < 22. Furthermore, our X-ray ERO 
samples arc similar in the proportions of starburst and AGN 
sources (on the basis of X-ray fluxes) and in their median 
hardness ratios of HR ~ 0-0.3. 

Vignali et al. (2003) repeated this cross-correlation on 
the same field, with the addition of more Chandra data 
bringing the exposure time to 2 Msec. Some additional 
EROs were detected, bringing the numbers to 10 of the 29 
HK' < 20.4 and 3 of the 9 if < 22 EROs. 



8.3 The X-ray-ERO Correlation at 9 > 

We find evidence, at the ~ 3.3a level, that Chandra sources 
are also cross-correlated with EROs at non-zero separations. 
The cross-correlation function, shown here (Figure 14a) for 



Figure 14. (a) Cross-correlation function, at non-zero separa- 
tions, between Chandra sources (73) and K < 21.5 EROs (179) 
on the ISAAC imaged area of the CDFS. This is also shown with 
the exclusion from the analysis of 17 Chandra sources identified 
as EROs. (b) Cross-correlation function between Chandra sources 
(73) and K < 21.5 EROs separated into pEROs (101) and ds- 
fEROs (78). 

Ks < 21.5 EROs, remains positive at 2 < S < 20 arcsec, 
where any signal would be the result of clustering between 
separate objects. A function 'A„(0~°-® — 12.24)' fitted at 
2.0 < ^ < 20 arcsec, gave an amplitude at one degree of 
40.8 ± 12.3 X 10"*, which is 68 ± 20 per cent of the ^(6*) 
amplitude (Table 2) for EROs at this limit. 

However, it seems possible that some of the signal at 
> 2 arcsec might be due to the fact that ~ i of Chandra 
sources arc also EROs, and EROs arc strongly clustered with 
each other. To investigate this, we excluded the 17 Chan- 
dra sources identified as EROs and cross-correlated only the 
remaining 56 with the ERO sample. The cross-correlation 
falls to zero at ^ < 2 arcsec (as expected), but it is little 
changed at ^ > 2 arcsec, where the best-fitting amplitude is 
52.4 ± 17.4 X 10"". This implies that, at « ~ 1-2, faint X-ray 
sources (AGN and starbursts) in general (i.e. not only the 
very reddened ones) cluster with EROs. 

We also calculated the cross-correlation of the Chandra 
sources with the separate pERO and dsfERO subsamples 
(Figure 16b), and fitted the same function as above at 2.0 < 
9 < 20 arcsec. At the Ks < 21.5 limit the best-fit A^, are 
34.5±27.0xl0"* (pEROs)and 39.2±24.2x 10"" (dsfEROs) 
- neither measurement is a 2a detection and there is no 
significant difference between the two, so a larger sample 
will be needed to investigate the dependence on spectral 
type. 

Almaini et al. (2003) cross-correlated bright sub-mm 
(SCUBA) sources and X-ray (Chandra) sources on the 
ELAIS:N2 field and detected a signal at 4o- significance, with 
a mean cross-correlation function of ~ 0.6 at non-zero sep- 
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arations of 5 < ^ < 100 arcsec. However, while the bright 
sub-mm sources are generally at high redshifts - e.g. Chap- 
man et al. (2003) find an interquartile range 1.9 < z < 2.8 
for a spectroscopic sample of 10 - the Chandra source N{z) 
is peaked at only z ~ 0.7 (Gilli ot al. 2003, from the CDFS 1 
MSec survey). Part of the cross-correlation could be due to 
gravitational lensing of the sub-mm sources by lower redshift 
large-scale structure, traced by the X-ray sources (Almaini 
et al. 2003b, and in prep.). 

On the other hand, the X-ray source N(z) also has an 
extended high-z tail, with ~ 20 per cent aX z > 1.5. Further- 
more these high-2 X-ray sources are clustered, with N{z') 
spikes at z — 1.618 and z = 2.572, and some of the faintest 
Chandra sources arc themselves SCUBA detections (Alexan- 
der et al. 2003). Hence, it is probable that a real clustering 
together of sub-mm and X-ray sources is also contributing. 

Fortunately, lensing effects will be much less significant 
in the case of our Chandra source - ERO correlation. Firstly, 
the effect of lensing on the observed surface density of any 
type of source depends on its number count slope; for a 
power-law ^ = 7 the cumulative count to magnitude m 
is modified by a muliplicative factor ix^'^'^~^ , where /x is 
the lensing magnification (Almaini et al. 2003b). The sub- 
mm number counts are extremely steep, 7 — 1, whereas the 
majority of the EROs are at Kg > 19.25 with 7 ~ 0.16. 
Lensing of ;u ~ 1.25 near a cluster would then enhance the 
probability of finding a sub-mm source by 40 per cent but 
reduce the rmmber of EROs by 10 per cent. Secondly, the 
redshifts of EROs (e.g. Cimatti et al. 2002a, 2002b) show 
much more overlap with the Chandra sources than do the 
sub-mm sources. Hence, we interpret our cross-correlation 
as evidence that these galaxies do trace the same large-scale 
structures. 



9 SUMMARY AND DISCUSSION 

1. We identify a sample of 198 EROs, defined here as 
7775 — Ks > 3.92 galaxies, to a limit Kg = 22 on public 
ESO/GOODS data covering 50.4 arcmin^ of the Chandra 
Deep Field South (CDFS). Of these, 179 are brighter than 
an estimated completeness limit of Kg = 21.5. 

The number counts of EROs flatten markedly at Ks ~ 
19.0-19.5, from 7 ~ 0.59 ± 0.11 to 7 ~ 0.16 ± 0.05. As we 
previously reported, counts of EROs are significantly lower 
than predicted by a model in which all present-day E/SO 
galaxies have undergone pure luminosity evolution (PLE). 
At Ks > 20, the ERO counts fall below even a non-evolving 
E/SO model, suggesting that the comoving number density 
(and not only the mass or luminosity) of red/passive galaxies 
is lower at z > 1. 

The ERO counts can be fitted much more closely by 
our 'merging and negative density evolution' (M-DE) model 
(Paper 1), in which a local luminosity function for E/SO 
galaxies is evolved through a combination of (i) passive L* 
evolution, (ii) merging at an evolving rate based on an ob- 
servational estimate (Patton et al. 2002), and (iii) a gradual 
decrease with redshift in the comoving number density of 
red/passsive galaxies, parameterized as R^. The best-fitting 
~ -0.49(±0.06) corresponds to a 42/61/68 per cent re- 
duction in (j)* at z = 1/2/3. This could be interpreted as a 
~ I fraction of the present day comoving number density 



of E/SO galaxies forming at « > 3 and the remainder from 
mergers and interactions of bluer galaxies (e.g. spirals) over 
all intermediate redshifts. 

2. We investigate the clustering of the CDFS EROs and 
detect a > 3cr signal in the angular correlation function, 
u){6), at Ks = 21 — 22 limits. The u!{9) amplitudes, com- 
bined with those from previous studies of ERO clustering 
(Paper I; Daddi ct al. 2000; Firth et al. 2002), arc inter- 
preted using models based on the Limber's formula integra- 
tion of our 'M-DE' model N{z) for EROs. The w(6») scal- 
ing of EROs at 7^ = 18-22 limits is well-fitted with either 
comoving (e — —1.2) clustering and a correlation radius 
ro = 12.5(±1.2) h^^ Mpc, or stable clustering (e = 0) and 
ro = 21.4(±2.0) h~^ Mpc. Our sample appears to be insuf- 
ficiently large for us to either confirm or exclude the recent 
claim that pEROs are more clustered at a ~ 1 than dsfEROs 
(Daddi et al. 2002). 

For full 7f -limited samples of galaxies, u!{6) is consistent 
the locally measured correlation radius of ro = 5.85 
Mpc, with little or no clustering evolution, e ~ —0.4-0. In 
contrast, EROs are even more clustered than local giant 
(L > L*) ellipticals (with ro — 8 h~^ Mpc), implying that if 
the EROs are direct progenitors of the E/SOs, they undergo 
strong clustering evolution (e < —1.2). 

This could be part of the same process as the 0* evolu- 
tion, if the oldest EROs form as strongly clustered, massive 
sub-mm sources, while the younger EROs, added at lower 
redshift, form from less strongly clustered disk galaxies and 
progessively dilute the red-galaxy ^(r) down to its present- 
day value. This scenario may be supported by a a prelimi- 
nary 'COMBO 17' finding (Phleps and Meisenheimer 2002) 
that E-Sb type galaxies show a rapid increase of clustering 
with redshift (e = -1.35 ± 0.24) over the < « < 1 range. 

The form of ERO clustering evolution should be bet- 
ter constrained in the near future as our planned spectro- 
scopic survey of these and other EROs will allow investigSr 
tion of ERO clustering in three-dimensions, and forthcoming 
clustering estimates for the sub-mm galaxies (Percival et al. 
2003) will provide an essential high-z point of comparison. 

3. There is an overdensity of EROs centered on the Chan- 
dra X-ray source XID:58, a Ks = 20.94 galaxy and itself 
an ERO, at R.A. 3''32'"11.85", Dec -27:46:29.14. Within 20 
arcsec of this object, we find 10 Ks < 22 EROs, compared 
to 1.37 expected by chance. Of these, 7 or 8 have colours 
consistent with z ~ 1.44, a recent photometric estimate for 
the Chandra source. 

The point-source X-ray emission from XID:58 corre- 
sponds to a luminosity Lx (0.5-10.0 keV) ~ lO'*^ '' ergs s"^ 
for h = 0.7, which indicates the presence of an obscured 
AGN. The host galaxy is of a modest optical luminosity - 
we estimate Mr ~ —22.2, which in our model would be near 
L* at this redshift. Willott et al. (2002) found an L ~ L* 
pERO in ELAIS:N2 (object N2:28) with a very similar X- 
ray to TT-band flux ratio. The host galaxy appears compact 
and irregular and we consider it more likely to be a dsfERO 
type, probably a recent post-merger, than a giant elliptical. 

4. We investigate the numbers of close (i.e. probably inter- 
acting) pairs of EROs, using the method of Woods, Falilman 
and Richer (1995) with the same pair selection criterion 
{P < 0.05). Prom the close pair counts we estimate the num- 
ber of interacting ERO companions per ERO (to Ks = 21.5) 
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as 0.062 ±0.019. This will a lower limit as it is probable that 
some companions will be missed, e.g. due to merging of the 
images. A very approximate correction for this revises our 
estimate up to ~ 0.073 ± 0.019. 

Patton ct al. (2002) estimated that normal, optically- 
selected galaxies at a; ~ 0.3 had 0.0321 ± 0.007 interacting 
companions per galaxy, this evolving as oc (1 + z)^-^ out 
to ~ 1. Hence, our result suggests that the pair fraction 
amongst EROs is higher than found in normal, optically- 
selected galaxies at z < 0.3, but consistent with those at 2: > 
1, even for our higher estimate. In contrast, radio-selected 
faint galaxies (primarily very strong starbursts) have a much 
higher interacting fraction, ~ 30 60 per cent (Windhorst ct 
al. 1995; Roche, Lowerithal and Koo 2002a), reflecting the 
importance of tidal triggering in producing the highest star 
formation rates. 

If the dsfEROs are undergoing disk-to-elliptical trans- 
formations, we might expect at least this type of ERO to 
show an elevated merger fraction. However, we previously 
(Paper 1) estimated, from the high dsfEROs/pERO ratio, 
that the dsfERO stage has a mean lifetime of ~ 1 Gyr, much 
longer than the timescales of either physical coalescence or 
merger-triggered starbursts (~ 0.1 Gyr). Furthermore, sim- 
ulations (Bekki and Shioya 2000) predict the most dust- 
reddened mergers will be those occurring with a retrograde- 
retrograde geometry of the galaxy rotation axes, and that 
these will starburst at a late stage of merging and lack ob- 
vious tidal tails. 

On this basis, the majority of the mergers showing 
dsfERO colours may be already be in post-burst, post- 
coalescence stages,and in Section 7 would have been counted 
as single galaxies rather than as pairs. If this is the case, 
many of the dsfEROs may be asymmetric or show other 
morphological features indicative of recent merging. We shall 
investigate the morphologies of these EROs in a forthcoming 
paper (Caputi et al., in prep). 

5. We find 17 coincidences between the 73 Chandra X-ray 
sources and the 198 EROs (Ks < 22 I - Ks > 3.75) in our 
data area, and from this estimate that 8.1 ± 2.1 per cent 
of the Kg < 22.0 EROs are X-ray sources above the Chan- 
dra 942 ksec detection limit. Of the 17 detected EROs, we 
classify 7 and pEROs and 10 as dsfEROs. 13 of the 17 are 
detected in the Hard band and have total X-ray fluxes suf- 
ficient to require the presence of AGN. These have a wide 
range of X-ray hardness ratios consistent with AGN obscur- 
ing columns Nh ^ 10^^-®"^^-^cm"^. The other four EROs, 
detected only in the Soft band (all dsfEROs), may be non- 
AGN starburst galaxies. 

The numbers and properties of the X-ray detected 
EROs are consistent with and very similar to the results of 
Alexander et al. (2002) for the CDFN. The results of VignaU 
et al. (2003) suggest that doubling the it Chandra exposure 
to 2 Msec would reveal a further ~ 10 per cent of the EROs 
to be faint X-ray sources. 

6. We also find a cross-correlation (at ~ 3.3c7 significance) 
between Chandra sources and EROs at non-zero (~ 2-20 
arcsec) separations. We estimate the cross-correlation has an 
amplitude 68 ±20 per cent that of ERO ij{6). Our sample is 
insufficient to determine the relative contributions of pEROs 
and dsfEROs. 

The cross-correlation is not significantly reduced when 



the 17 Chandra sources coincident with EROs are excluded. 
This implies that EROs are clustered with faint X-ray 
sources at z ~ 1-2, ie. with AGN and powerful starbursts, in 
general and not only with those that are heavily reddened. 
As discussed in the Introduction, this is as expected for the 
scenario (Archibald et al. 2002; Granato et al. 2003) where 
sub-mm sources evolve into EROs via a QSO phase. 

We plan to use ultrardeep spectroscopy of these EROs 
and others to investigate further (and in three dimensions), 
the environments and cross-correlations of pEROs, dsfEROs 
and faint X-ray sources. 
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